Spin-transfer torque in magnetic tunnel junctions: Scattering theory 
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We study the bias-dependent spin-transfer torque in magnetic tunnel junctions in the Stoner 
model by scattering theory. We show that the in-plane (Slonczewski type) torque vanishes and 
subsequently reverses its direction when the bias voltage becomes larger or the barrier wider than 
material and device-dependent critical values. We are able to reproduce the magnitude and the 
bias dependence of measured in-plane and out-of-plane torques using realistic parameters. The 
condition for the vanishing torque is summarized by a phase diagram depending on the applied bias 
and barrier width, which is explained in terms of an interface spin polarization and the electron 
focusing by the barrier. Quantum size effects in the spin-transfer torque are predicted as a function 
of the thickness of a normal metal layer inserted between the ferromagnet and tunnel barrier. 



I. INTRODUCTION 

Magnetic tunnel junctions (MTJ) are layered struc- 
tures in which an insulating tunnel barrier (I) sepa- 
rates two ferromagnetic layers (F).^'^ The interplay be- 
tween electronic currents and an order parameter differ- 
ence, i.e. magnetizations rotated away from the equi- 
librium configurations, is the magnetic equivalent to 
the Josephson effect in superconductivity. The MTJ 
with a thin normal metal insertion layer is the only 
magneto-electronic structure that shows quantum size 
effects on electron transport."^ MTJs based on epitax- 
ial MgO barriers^ are used in the magnetic random- 
access memory (MRAM) devices that are operated by 
the spin-transfer ioxqne,^ MTJs have been studied 
vigorously, initially focusing on the tunnel magneto- 
resistance (TMR) ^i^iSiiS More recently, the focus shifted 
to spin-transfer torque and current-induced magnetiza- 
tion switching^aais^isaeiiasiQ^^^^^ On the 

theoretical front, spin-transfer effects have been studied 
extensively in metallic spin valve structures based on var- 
ious models i^^i^^'^^i^^'^^i'^° For tunneling structures, such 
studies are still relatively scarce J^ii^ii^iii^i^ 

Here we report a model study of the spin-transfer 
torque in magnetic tunnel junctions. Since the ferromag- 
nets are separated by tunnel barriers, we cannot use theo- 
ries existing for metallic structures that are mostly based 
on semiclassical methodsi^i^ Instead, we chose a fully 
quantum mechanical treatment of transport through the 
tunnel barrier by scattering theory. The high quality 
of MgO tunnel junctions and the prominence of quan- 
tum oscillations observed in FNIF structures (even for 
alumina barriers) provide the motivation to concentrate 
on ballistic structures in which the transverse Bloch vec- 
tor is conserved during transport. We qualitati vely (a nd 
even quantitatively) confirm the results in Ref. Il5lll9ll20l . 
However, our model is simpler and physically more trans- 
parent than the tight-binding method used in Ref. [l^ 
and the numerical studies of Refs. [l9| and Hfli We are 
able to reproduce simultaneously both the in-plane and 



out-of-plane torque experimental data using realistic ma- 
terial and device parameters, in contrast to a fit based on 
the tight-binding model. We also show finite zero-bias 
out-of-plane torque for asymmetric structures. Scatter- 
ing theory enables us to distill a clear physical picture 
of the peculiarities of the spin-transfer torque in MTJs, 
which allows us to understand why and when the torque 
goes to zero. The torque zero-crossing condition can be 
summarized by a phase diagram spanned by the applied 
bias and barrier width parameters. With our approach 
we can go beyond the FIF MTJ and study the effect of 
a normal metal insertion (FINF structures). Quantum 
size oscillations in the torque are predicted as a function 
of the thickness of the N insertion layer. 

This paper is organized as follows: Section [IT] intro- 
duces the FIF and FNIF structures and the scattering 
theory. In Section [III[ approximations are introduced in 
order to derive analytic expressions. Section [TVl presents 
our main results. Section |V] compares our model with 
experimental results. Section IVll and IVTIl contain a brief 
discussion and summary, respectively. Two appendices 
are attached at the end. 



II. STRUCTURE & METHOD 

We consider multilayers as shown in Fig. [IJa) in which 
two semi- infinite F leads (F(L) and F(R)) are connected 
by an insulating layer (I) of width d and a non-magnetic 
metal layer (N) of width a. The magnetization direction 
of F(L)/F(R), mi/m2 (|mi| = |m2| = 1), is treated as 
fixed/free. This structure reduces to a conventional FIF 
MTJ when a = 0. 

Let A, B, C, D, C , D', E, F be the spin-dependent am- 
plitudes (At = (A|, A|)) of flux- normalized spinor wave- 
functions at specific points. The scattering states can be 
expressed in terms of two incoming waves A and F, such 
as: 

C = .sc'aA + sc'fF, (1) 




rent reads, 



FIG. 1: (Color online) (a): FNIF heterostructure, in which 
Si^2 indicate two different interface scattering regions; (b): 
The potential profiles (at positive bias) for majority and mi- 
nority electron spins in F are shown by solid and dashed lines 
respectively. The exchange splitting is A and the tunnel bar- 
rier has height Uh relative to the Fermi energy Ep. The ap- 
plied bias V pictured in (b) corresponds to a net electron 
(particle) flow from right to left. 



where sq'A and sc'f are 2x2 matrices in spin space 
that can be constructed by concatenating the scattering 
matrices of region 5*1 2 and of the insulating layer bulk 
(see Fig. [1]). To leading order of the transmission (ib) 
through the bulk I (similar expansions hold for sd'A and 
sd'f as well): 



q,q' 



(3) 



4e, 



-Tr, 



Im (^seaSea/l - sefSef/r 



where Tr^ [• • •] denotes the spin trace, and the sum- 
mation is over all the transverse modes at energy E. 
/l = fhiE) and /r — fB.{E + eV) are (zero tempera- 
ture) Fermi-Dirac electron distribution functions in the 
left and right reservoirs. We are therefore disregarding 
any spin-accumulation in the ferromagnet, which is vahd 
for tunnel junctions of current interest in which the spin- 
flip rate in the ferromagnet is larger than the tunnel rate. 
The scattering matrices depend on V by the bias-induced 
potential profile. The spin current, or the angular mo- 
mentum current, at the left side of the I/F(R) interface 
(within I) reads 
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(4) 



Since the spin current deep in the FM lead is longitudi- 
nal to the magnetization, the torque N acting on F(R) is 
equal to the transverse component of the incoming spin 
current that is absorbed at the interfacei^ i^^'^^i^^ 



N = Js - (Js • m2)m2 = N|| 



Nj 



(5) 
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where ii, 2/^^1,2 are the 2x2 transmission/reflection ma- 
trices for Si^2 (see Fig. [T|) and tb/fb are the spin- 
independent transmission/reflection coefficient for the in- 
sulating bulk material that are proportional to the unit 
matrix in spin space and therefore without hat. The 
primed and unprimed versions indicate scattering of elec- 
trons impinging from the left and right, respectively. The 
reflection coefficient rb is due to the impurity scattering 
inside the bulk insulator, which, as ib, contains an expo- 
nential decay factor representing evanescent states in I. 
For this reason the magnitude of rb is comparable to or 
much smaller than that of ib depending on the density of 
the impurities. All scattering matrices are matrices in k- 
space defined by the propagating states of left and right 
leads in the energy window available for transport, la- 
beled by their transverse wave vectors in the leads: q, q' 
(the band index is suppressed) at a given energy. 

An applied bias voltage V drives a (conserved) charge 
current Jc and (spatially-dependent) spin current Jg 
through the device. At zero temperature, the charge cur- 



with the in-plane (Slonczewski) torque N|| cx m2 x (mi x 
m2) and out-of-plane (field-like) torque oc rtii x m2. 
Similarly n = jg - (js • m2)m2 = n|| -I- n_L. 

At low bias, the non-equilibrium part of the spin cur- 
rent is proportional to the bias voltage Jg — Jg = GsV, 
where Jg is the equilibrium spin current that is related 
to interlayer exchange coupling at equilibrium, and Gg is 
the spin conductance: 
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where the scattering matrices are evaluated at zero bias 
V — 0, and the summation is over transport chan- 
nels at the Fermi energy. We define the linear-response 
torkance^" T = N/V, and T = Gs - (Gs • m2)m2 = 
Til + T_L and r = gs - (gs • m2)m2 = rii + t±. 



III. APPROXIMATIONS 

We assume in the following that the spin is conserved 
during the scattering. Then ti for Si {i — 1, 2, similar for 
fi) is diagonal when choosing rrii as spin quantization 
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axis: Expanded in Pauli matrices a = (o-^;, CTj,, iTz), ti = 
tf + t^a ■ m„ with tf = (4 ± t\)/2. if {a =T, i) is the 
transmission amphtude for spin-a for spin quantization 
axes m.i through the scattering region Si. In the absence 
of impurities (rb = 0) and to leading order of tb^ 

sc'A^th{tt +ti^ -mi), (7a) 
SD'A = {rttbfl + rj^tb^r™! • ma) (7b) 
+ (T • {r^^t^xni + r^ibtr™2 - ir^tbt^nni x m2), 

Similar expansions hold for sc'F and sd'F- 

Next, we adopt the free electron approximation tai- 
lored for transition metal ferromagnetsi^ We assume 
spherical Fermi surfaces for spin-up and spin-down elec- 
trons (in both F(L) and F(R)) with Fermi wave-vectors 
kl = y/2rriE^^ and = ^j2m{EY - A)/h'^, with an 
effective electron mass m in F. Electrons in N are as- 
sumed to be ideally matched with the majority electrons 
in F {kp = fcp). The effective electron mass in the tun- 
nel barrier is assumed to be Wb = (3m. Fig. [TJb) shows 
the adopted potential profile with barrier height t/b • We 
assume an applied potential bias that is smaller than the 
barrier height \eV\ < U\y that fully drops over the tunnel 
barrier. Positive bias corresponds to charge current flow 
from left to right and electron particle flow from right to 
left. As in Refs. 15,18,19, we assume that energy and 
transverse wave-vector q are conserved, thus disregard 
any impurity/interface roughness scattering. All scatter- 
ing matrices then become diagonal in k-space. 

In the free electron model, the flux-normalized wave- 
functions in N and F are: 



in N and 



/ho 

A.™ 



in F, (8) 



Using the WKB approximation,^^ the wave- function in 
the tunneling barrier is 



^ ~ 'U I 



(9) 



with 



^(u^ + EF-E-eV^)+q^, (10) 



where E is the energy of the electron. The WKB ap- 
proximation is valid when the potential profile varies 
slowly in space within the tunneling barrier, i.e. k'{x) <C 
K^{x). The transmission coefficient through I reads 

tb — exp[— K{'w)dw\. 

For finite bias, from Eq. (g]) and Eq. 
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= 2tl Re (r-r-fi^ + T-r[ /r) rm x ma, (lib) 
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where = \t]\'^ + \tl\'^ is the average transmission prob- 
ability for scattering region Si, T~ = p^T^ = \tj\'^ - |tf p 
with polarization pi = T~ /T^, and /o is the equilib- 
rium distribution function at zero bias, in Eq. (jllb|l 
includes both the equilibrium and non-equilibrium con- 
tribution to the out-of-plane torque. The former {n\ 
in Eq. (jllc|) ) is related with the non-local intcrlaycr 
exchange coupling. '^^ The non-equilibrium contribution 
is therefore nY = — n*|. The optical theorem 
2Im (^'j^a) ~ "^1 2 ('^'^e Appendix A) is used in the deriva- 
tion of Eq. (jllap to get rid of all internal reflection in I. In 
Eq. (11a), we observe that the in-plane torque is caused 
by the polarization of the current at the left interface 
that is expressed by . The subsequent absorption of 
the spin current by the second magnet is governed by the 
geometrical projection expressed by the vector product 
and the total transparency of the second interface T^. It 
follows from Eq. (lib) that the out out-of-plane torque 
has a very different origin. It does not depend directly on 
the difference of the electron distributions on both sides 
of the junctions, but consists of two independent contri- 
butions from both reservoirs. Each contribution consists 
of the spin-polarization of the first interface, but is sensi- 
tive to the phase of the reflection coefficient of the second 
interface. The out-of-plane torque can be interpreted as 
the net spin created at one interface that while reflected 
at the second interface briefly precesses in the exchange 
field of the second ferromagnet. 

With vanishing bias, ib = exp(— kc?) with k = 
^2mbC/b/?i2 + g2. By Eq. jS]) and Eq. 
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T2+m2 X (mi x ma), 



(12) 



and T ^ = 0. For reference, the conductance within the 
same theoretical framework is given byi^^ 



2h 



(13) 



The vector product \xn.2 x (mi x ma) | = sin 9 in Eq. (fl3 
and mi • ma — cos 9 in Eq. leading to the well- 

known geometrical dependence of the angular transport 
properties of tunnel junctions.^ The vanishing of the out- 
of-plane torque, T± = 0, is a rather general result that 
holds for symmetric tunneling junctions and spin valves 
in the linear-response regime We consider a symmetric 
system with an applied voltage —V/2 to the left and a 
voltage V/2 to the right reservoir. To the second order in 
the bias voltage, the spin-current in the spacer between 
the ferromagnets can be expanded as 

Is = [Aim.1 + Bima + Cimi x ma] V 

+ [Aami +B2ma + C2mi X m2]F^ (14) 

When applying the mirror operation left ^ right (1 ^ 
2,— V/2 ^ V/2, Is ^ —Is) symmetry requires that 

-Is = [Aim2 + Biuii + Cima x mi] {~V) 

+ [Aama + B2mi + C2m2 x mi] {-Vf, (15) 
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which should be identical to Eq. (|T4|) . 
Bi,Ci = 0,^2 = —B2, whereas C2 
Then the torque on m2 is 



Therefore Ai = 
is not restricted. 



in-plane 



out-of-plane 



N = m2 X (I, X m2) (16) 
= [AiV - A2V'^] m2 X (mi x m2) + C2V^^mi x m2. 

This proves that, for symmetric systems, the out-of- 
plane torkance (cx 2C2V) vanishes at = 0. It also 
shows that beyond linear response, there are quadratic 
(in bias) contributions to both the in-plane and out-of- 
plane torque. The argument does not hold for asymmet- 
ric tunneling systems. An experimental zero-bias out-of- 
plane torkance should therefore provide interesting infor- 
mation on MTJ asymmetries. 



IV. RESULTS 

In this section, we discuss three different structures: 
(A) a symmetric FIF magnetic tunneling junction, (B) 
an asymmetric FIF structure in which the left and right 
FM layers have different exchange splitting, and (C) an 
FNIF structure in which a non-magnetic layer is inserted 
between the insulator barrier layer and one of the ferro- 
magnetic layers. 



A. Symmetric FIF 



For a symmetric Fe/MgO/Fe MTJ: kl = 1.09 A"^ and 

4 0.42 for Fe,-^ and ^ 1.0 - 1.2 eV and /? 
rrib/m = OAm for MgO^i^ii^i^e This implies Ep ~ 4.5 
eV, A ~ 3.8 eV « 0.85£;f, and Ub « 0.25£;f. For an 
FIF structure (a — 0), both Si and 5*2 contain only a 
single interface. Using the potential profile in Fig. [Ijb), 
we have 
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kl + iK{Q)/l3' 
~k1 + iK{Q)/l3 
k'[ + iK{0)/(3 ' 



H — 



2y/ik-K(d)/(3 

k^ + in{d)/(3 ' 

-k^ +iK{d)/f3 

k^ +iK{d)/l3 ' 



(17a) 
(17b) 



where a =t, i and fc^+gr^ = 2mE/fi^, k\^+q^ = 2m{E- 
A)/h'^. fcf are defined similarly with E replaced by _E + 
eV. We set t^ ^ = when Im (kl ^) ^ 0. 

Fig. [21 shows the computed bias dependence of the 
in-plane torque N|| (left) and the non-equilibrium part 
{i.e. not containing the equilibrium interlayer exchange 
coupling) of the out-of-plane field-like torque N'^(l/) = 
Nj^(F) — N^(0) (right) at various exchange splittings 
for TOb — 0.4m at d = 1.0 nm. The equilibrium exchange 
coupling gives rise to an effective magnetic field in the 
LLG equations that we do not explicitly discuss. The 
main features of these curves are: 1) the in-plane torque 
has both linear and parabolic contributions; 2) the field- 
like torque is parabolic like. These plots are very similar 
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FIG. 2: (Color online) The magnitude of Ny (left) and N'^ 
(right) acting on the right magnetization in Fig.[T]vs. applied 
bias eV for (3 = 0.4 and d = 1.0 nm. Inset figure shows the 
out-of-plane torque including the equilibrium contribution at 
zero bias. 



to the corresponding plots by Theodonis et alf^ mean- 
ing that the band structure effects caused by the tight- 
binding approximation are not important. The in-plane 
torque at negative bias and small positive bias is "nor- 
mal" , but changes sign at higher positive bias, where nor- 
mal means that the direction of the torque in FIF is the 
same as the torque in metallic spin valves predicted,^ i.e. 
the torque curve appears in the second and fourth quad- 
rant in the left panel of Fig. [2l When the torque curve 
is found in the first or third quadrant, we say that the 
torque is reversed. In the normal region, the positive bias 
(electron flow from right to left) favors the anti-parallel 
configuration and a negative one the parallel configura- 
tion. In the reversed region, on the other hand, a current 
polarity that stabilizes the parallel configuration in the 
normal region has the opposite effect. 

The zero-crossing of the in-plane torque in Fig. [2] can 
be traced to the sign change of in Eq. piap . i.e. 
the sign change of the polarization of Si (the F(L)/I 
interface)pi = /T^ < OM- The polarization pi cx 

(0) — k-^ki, which can take any sign depending on pa- 
rameters choseni^ (see Appendix B for a more detailed 
discussion of this point). The vanishing torque phe- 
nomenon becomes more transparent without an effective 
mass mismatch, i.e. for (3=1 instead of /? = 0.4 used in 
Fig. O The polarization vanishes when 



K^iO)^k^k^=0. 



(18) 



Since k(x = 0) (near left interface) increases and k^ de- 
creases with q, Eq. is fulfilled at a certain critical 
value gc- The latter increases with the electron energy 
E, because k(0) decreases and increases with E. This 
can be seen clearly from the following equation: 



2m 



E- 



{Ef + £/b)^ 
2{Ef + Ub) - 



(19) 



4 



(a) 



® 




rp+ 




t\eV\ 




■i-- 

1— 


-- 

— f— 







(b) 










+ 


+ 











. V 


T+ 

\ Y 


^ (S> 








-i--- 

— f 


- -i - 
— f — 





(7b 



negative bias: V < positive bias: V > 

FIG. 3: (Color online) Positive and negative bias situation. 



which impUcs that at the Fermi energy qc is well defined 
(g^ > 0) only when < Ef{Ef - A). 

When < at low bias there is no polarization sign 
change for any q, and the (in-plane) torque behaves nor- 
mally. When the potential profile becomes distorted by 
an applied bias as in Fig.[3Jb), the electrons injected from 
the right lead have a maximum energy E = Ep + eV. 
When the applied bias V is large enough, we reach the 
regime q^ > 0, and a polarization sign change of the left 
interface comes into play. As increases further, more 
and more electron contribute to the opposite torque. 
When V is large enough, the total torque changes sign as 
seen in Fig. [2l On the other hand, when the applied bias 
is negative (see Fig.|3Ja)), the transport is dominated by 
the electron injected from the left lead. The electron en- 
ergy and effective barrier height at the left interface do 
not change with applied bias, which means Tj" (and so 
the polarization) does not change either. Therefore, we 
do not see a zero-crossing (on the right magnetization) 
at negative bias in Fig. [2l 

We can analogously understand the dependence of the 
in-plane torque on the barrier width d. We know that 
there is a polarization sign change for q < qc if Ub is not 
too high (such that q^ > 0). In timneling junction, the 
transport is dominated by electrons with small q because 
of the focusing effect due to the exponential extinction 




-0.1 0.1 0.2 
bias voltage: eV/Ep 

FIG. 4: (Color online) Critical barrier width dc vs. barrier 
height Uh at linear-response (left) and applied bias V (right). 
A = 0.85£;f. Critical width for TMR at /? = 0.4 (thin solid 
black line with "-|-" symbols) is also shown in the right panel. 



factor in th- When the electron with q values smaller 
than qc dominate, the torque reverses sign. At a critical 
barrier width dc, the contributions from q < qc cancel 
those from q > q^, and the torque or torkance vanishes: 
N||((ic) = T||((ic) = whereas it has opposite direction 
for d > dc- The left panel of Fig. 2] shows dc vs. Uh in 
linear-response for various values of the barrier effective 
mass, dc increases with Uh simply because the polariza- 
tion sign change behaviour is less prominent at higher 
barrier heights (q^ decreases with Uh)- 

The right panel of Fig. 0] shows the bias dependence 
of dc at barrier height Uh = 0.25£'f for various barrier 
effective masses. Because both the positive and nega- 
tive bias reduces the average barrier height, th increases 
with \V\. However, the applied bias is very ineffective 
in changing the focusing behaviour, i.e. the transmis- 
sion is hardly less focused by the reduced average barrier 
height. This is very unlike the geometrical barrier width 
d, to which the focusing very sensitively: larger (smaller) 
d means more (less) focusing. Rather, a positive bias 
enhances the polarization sign change behaviour, which 
leads to a smaller critical barrier width dc. Since there is 
no polarization sign change (Qc < 0) at zero or negative 
bias for (3 — 0.7 and 0.4, a torque zero-crossing is not 
observed. For f3 = 1.0, polarization could change sign 
at zero bias, hence the torque zero-crossing is also ob- 
served at negative bias. As mentioned before, negative 
bias does not change or polarization of the left inter- 
face, but it does change . At negative bias, the barrier 
height at the right interface is reduced by \eV\, which 
leads to the decrease of at small q, thus the polar- 
ization sign change behaviour is weakened because of the 
smaller product T^T2 at small q where the product is 
negative. Weaker polarization sign change then requires 
larger dc for torque zero-crossing at negative bias and 
we see dc increasing with negative bias for [3 — 1.0. For 
comparison, a critical barrier width for the sign change of 
TMR is also calculated for (3 — 0.4, and is shown as the 
solid black curve with "+" symbol in the right panel of 
Fig. m Since the TMR is symmetric in the applied bias 
for the symmetric structures, the curve is also symmetric. 



B. Asymmetric FIF 

As discussed at the end of Sec. IIII[ the zero-bias 
(linear-response) out-of-plane torkance does not vanish 
for asymmetric structures. Fig.l^lshows both the in-plane 
and out-of-plane torkance for the right FM layer at zero 
bias [V — 0) for an asymmetric FIF structure, in which 
the left and the right FM layers have different exchange 
splitting: Ai — O.SSi^F for the left and A2 for the right, 
where the latter varies from to E-p. From Fig.O we can 
see that the out-of-plane torkance is generally non-zero 
for asymmetric structures when A2 9^ Ai, and it vanishes 
when the right layer becomes non-magnetic (A2 = 0) or 
when the structure becomes symmetric (A2 = Ai). The 
in-plane torkance in Fig. [5] decreases with A2 simply be- 
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FIG. 5: Linear-response torkance vs. exchange splitting A2 
in the right FM (Ai = O.SSSf, C/b = 0.25£F,d = 1.0 nm). 



cause the average transmission probability through the 
right interface decreases. 



C. FNIF 

An FNIF structure, with a non-magnetic layer of width 
a between one of the F layers and the insulating I layer, 
has never been studied in the regime of spin-transfer 
torque. Such an asymmetric FNIF device can be op- 
erated in two non-equivalent modes (as in any non- 
symmetric MTJ): the left F is static and the right F (F) 
is free (mode 1: FNIF), and vice versa (mode 2: FNIF). 
Eqs. pT|) - apply to mode 1, and apply to mode 
2 with subscript 1 and 2 exchanged. The a-dependence 
of the in-plane torkance (in linear-response) is shown in 
Fig. [S] The sign of the in-planc torkance can be con- 
trolled by a in mode 1, but not in mode 2. This sign is 
determined by the sign of . In mode 1, T-^{a) covers 
region F(L)-N-I and its sign can be modulated by the N 
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FIG. 6: (Color online) Linear-response in-plane torkance vs. 
non-magnetic insertion layer width a (A = 0.85Ef ,Uh ~ 
0.25EF,d = 1.0 nm). 



insertion layer width a. However, in mode 2, Tj~ covers 
I-F(R), which is independent of a, therefore the sign is 
unchanged. The a-dependence of the in-plane torkance 
in mode 2 comes from T^{a), which is always positive. 
Due to the aliasing effect caused by discrete thickness of 
the N layer,^ the period of the quantum oscillation in 
Fig. [5] should be about Tr/\k-p ~ n/ \\ instead of tt/Zcf ~ 3 
A shown in the figure, where A is the monolayer thickness 
for N layer. 

The asymmetry in FNIF structure shall also give rise 
to a finite linear-response out-of-plane torkance, which 
shows similar oscillations as the in-plane torkance in 
Fig. [SI The magnitude of the zero-bias out-of-plane 
torkance could be comparable to the in-plane counter- 
part. 



V. COMPARISON WITH EXPERIMENTS 

Using realistic material parameters and the geometry 
parameters provided by Ref. we are able to repro- 
duce, even the absolute scale, the experimental data from 
Ref. [l^ as shown in Fig. [71 which includes the bias de- 
pendence of the in-plane torque N||(y) and the non- 
equilibrium part of the out-of-plane torque N'^(y) — 
Nx(y) — Nj^(O) and the corresponding torkance. The 
experimental data for the torkance in the bottom pan- 
els of Fig. [7] are adapted from Fig. S3(d) (/^stft) ^^'^ 
Fig. 2 (I-V data) of Ref. HI by {dN^^ ^/dV)/ sin 9 = 
{dI/dV){dNn^j^/dI)/ sine = (d//dl/)/3^^ p^. Our model 
appears to have a problem with the upturn of the torque 
at higher positive bias. The fit in Ref. [12 based on 
the tight-binding model of Ref. tlSi is slightly better in 
this respect for a rather large exchange splitting. How- 



in-plane 



out-of-plane 




-0.2 0.2 0.4-0.4 -0.2 0.2 0.4 
bias voltage: V (Volts) bias voltage: V (Volts) 

FIG. 7: (Color online) Fitting (solid curve) of the experi- 
mental data (dots) from Ref. [23. Top: in-plane torque 
AA'^II (left) and out-of-plane torque ANi_ (right); Bottom: 
in-plane torkance A{dN^^/dV)/ sin 9 (left) and out-of-plane 
torkance A{dN'j^/dV) / sinO (right). The following parame- 
ters are used in all fittings: Ef = 4.5 eV, A = O.STBf ^ 3.9 
eV, Uh = 0.23-Bf ~ 1.0 eV, P = 0.36; d = 1 nm, cross section 
area A = 70 nmx250 nm, and 9 = 137° from Ref. lH. 
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ever, the out-plane torque is poorly reproduced for the 
same parameter set. In contrast, we succeed with a single 
set of (realistic) parameters to reproduce both in-plane 
and out-of-plane torques (torkance). The resistance for 
this particular device in our model is R{9 — 137°) « 
150 r2, which is consistent with the experimental values 
(~ 200 n). However, the TMR value in our calculation 
(about 15%) is considerably smaller than that in the ex- 
periment. We believe that the spin-transfer torque is 
better represented by the free electron model than the 
TMR because TMR cx piP2 whereas torque on m2 oc pi . 
If the polarization pi , p2 at the interfaces are underesti- 
mated by a factor of rj, the TMR value is too small by a 
factor of 77^, whereas the spin-transfer torque is affected 
only by a factor of 77 (see also Section |VT|) . In addition to 
this, the TMR value depends sensitively on the exchange 
splitting A. For instance, the TMR value increases from 
15% to 30% when A increases from 0.87£;f to 0.9£;f, 

Another set of the experimental data is shown in Fig. [8] 
adopted from Fig. 3(a) of Ref. H!]. The experimental 
data are now the in-plane (red squares) and out-of-plane 
(blue diamonds) torkance. We fit the out-of-plane and 
the linear part of the in-plane torkance again on absolute 
scale. Theory contains a quadratic component of the in- 
plane torque, which does not show up in this experiment. 
The resistance for this particular device in our model is 
R{9 = 71°) « 4 kil consistent with the experimental 
value (^ 3.5 kil), whereas our TMR 5% again is too 
small. Note that the torkances are much smaller than in 
Ref. 22, because of the thicker barrier and the smaller 
cross section area, as reflected by the higher resistance. 

Deac et al^ also measured the in-plane and out-of- 
plane torque in a MgO based tunneling junction. The 
out-of-plane torque in this experiment agrees well with 
other experiments and theory. However, the in-plane 




bias voltage: V (Volts) 

FIG. 8: (Color online) Fitting (curves) of the experimental 
data (dots) from Ref. I21I . The following parameters are used 
in both fittings: Ef = 4.5 eV, A = 0.85Ef « 3.8 eV, Ub = 
0.25£'p ~ 1.1 eV, /3 — 0.43; d — 1.25 nm, cross section area 
A = 50 nmxlOO nm, and 6 = 71° from Ref. |2l| (Notice the 
sign convention in Ref. [U is opposite to that in Ref. [2^ . 



torque depends parabolically on the applied bias {AV^ 
for positive bias and —BV^ for negative bias, where A 
and B are positive constants), which is quite different 
from both Ref. [2l| and ll^. A voltage noise measure- 
ments done by Petit et al^ also suggest a linear out-of- 
plane torkance (or parabolic out-of-plane torque), which 
is about 20% of the in-plane counterpart. Hence all 
experiments and theories appear to agree on the out- 
of-plane torque, whereas consensus about the in-plane 
torque has not been reached yet. 



VI. DISCUSSION 

Because of the high quality of epitaxial MgO tunnel 
layers )^ we ignored interface roughness and barrier dis- 
order. The main effect of the geometric interface rough- 
ness is to reduce the nominal thickness of the barrier—. 
Impurity states in the barrier generally increases tunnel- 
ing because of the opening of additional tunneling chan- 
nels with lower barrier height U!^ <Uh- Impurities states 
also weaken the spin-dependent effects when spin-flip is 
involved. In general, interface roughness and disorder 
can be important quantitatively, but have been shown 
not to qualitatively change the features predicted by a 
ballistic model^a^^i^ 

The free electron Stoner model is only a poor represen- 
tation of the real electronic structure of transition metals 
for the tunneling problem: it fails to properly reproduce 
the nearly half-metallic nature of transition metal ferro- 
magnets based MgO tunnel junctions, that is caused by 
the symmetry of the bands at the Fermi energy ,~ leading 
to the underestimated TMR ratios by our model noted 
above. On the other hand, the band structure calcula- 
tions by Heiliger et alr^- show that the free electron model 
can perform quite well as far as the torque is concerned. 
We explained this in Section |TV] by its dependence on 
only one interface polarization leading to a better per- 
formance of a model that is not accurate in this respect. 
The band structures calculations in Ref. [l3 indicate that 
the torque is strongly localized to a few monolayers which 
is in support of our simple model. 

The issue of the wave-function symmetry should also 
be considered when a normal metal is inserted. When 
the electrons with wave vector normal to the interface 
dominate, the normal metal (Cr) is actually a potential 
barrier for the Fe majority spins*'^ rather than a potential 
well as assumed here. 

In contrast to metallic spin valve structures, in which 
the out-of-plane torque is generally less than 10% of the 
in-plane counterpart, the out-of-plane contribution has 
been found quite large in tunneling junctions (a 30% con- 
tribution at high bias is measured in Ref. 1211 ). Close to 
the zero crossing of the in-plane torque at positive bias 
the out-of-pane torque should become dominant. 

An experimental "phase diagram" that can be com- 
pared with Fig. 2] would constitute a stringent test of our 
predictions. Since the barrier height and the effective 
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electron mass in the barrier can not be controlled, we 
suggest to measure the torques systematically for several 
MTJ structures with different barrier width (otherwise 
identical) to test the red curve in the right panel of Fig. H) 
The zero crossing of the in-plane torque is predicted to 
occur at voltages that are too high for the experiments 
m Refs. [lUilil. For wider tunneling barriers it should 
happen at smaller voltages. 



VII. SUMMARY 

To summarize, scattering theory of transport is used 
to calculate the spin-transfer torque for a magnetic mul- 
tilayer structure at finite bias. The experimental spin- 
torque data (both in-plane and out-of-plane) can be re- 
produced using realistic parameters in our model. The 
spin-transfer torque on a given layer may change sign for 
only one bias polarity. The torque zero-crossing is caused 
by the combined effects of the polarization sign change at 
the FI interface and the focusing effect of the barrier. The 
bias voltage required for torque zero-crossing decreases 
as a function of the barrier width. The out-of-plane 
torkance at zero bias vanishes for symmetric FIF struc- 
tures, but remains finite for asymmetric structures. In 
FNIF structure we find on top of the previously reported 
oscillating TMIJiSi and charge pumping voltageS that the 
spin-transfer torque also oscillates and may change sign 
with the N layer thickness. 



waves e^^^^ / Vkn in N with mode index n, and exponen- 
tial solutions e^'*'"^/y'K^ in I with mode index m. The 
electron with wave function e^'^^^/y^K^ in I is reflected 
with amplitude Vmp^ and transmitted into N with ampli- 



tude tnp. Then the wave functions in I and N are: 



np 



(Al) 



The current in I is given by the imaginary part of the 
reflection amplitudes since 



2h 



Ini(rpp) , 



m ' ^ 



The current in N reads 

/n = - 5] Im («*9.«) = - ^ |i„p| 
m ^ — ' m ^ — ' 

n n 

By current conservation: /j = In, we have 

n n 

This relation reduces to 

2 Im (t^pp) — Kppl = Tpp. 
for the ballistic model used in the text. 



(A2) 



(A3) 



(A4) 



(A5) 



APPENDIX B: POLARIZATION SIGN CHANGE 
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APPENDIX A: OPTICAL THEOREM IN 
TUNNEL JUNCTION 

The scattering matrices in Eq. ^ include the reflection 
amplitude inside the insulator (f'l and f2). Expressions 
become more transparent when the reflection amplitudes 
are replaced by the transmission probabilities, however. 
This can be achieved by the optical theorem tailored for 
tunnel junction that reflects current conservation. Note 
that the equivalent statement in metallic systems is the 
well-known relation |rp -I- |ip — 1 that follows from the 
unitary of the scattering matrix. 

The optical theorem for light is derived from conserva- 
tion of energy, whereas in electronic transport it is based 
on conservation of charge. We consider here the non- 
standard situation of the interface between a metal and 
a tunneling barrier, for which the unitary of the scatter- 
ing matrix cannot be invoked without some care. Let us 
consider a non-magnetic IXN structure, where X could 
be basically anything. We assume flux-normalized plane 



For a better understanding of the polarization sign 
change, let us inspect the simple potential barrier de- 
picted in Fig. O (thick solid Hue), ignoring the spin de- 
pendence for the moment. The barrier width is d and 
the relative barrier height Uh — h^K^ /2m. As seen in 
Fig. [9l Uh is measured relative to the longitudinal elec- 
tron energy in the barrier E — E± = E^^ + V, where E is 
the total electron energy and -Ey and E± are the longi- 
tudinal (normal to the interfaces) and transverse kinetic 
energy, and V is the band edge. By solving this standard 
quantum mechanical exercise, we find the transmission 
probability through the barrier 



T 



l+(4L±i^sinh2(^d) 



{E\\ + Uhf 



(Bl) 



z 




^2 
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FIG. 9: (Color online) Potential barrier. 
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0.01 








FIG. 10: T vs. E^JUi ((7b = 0.25£f,-B = Ef,I3 = 0.4, and 
d — 1 nm). 

where the approximation is accurate when nd ^ 1. 
Eq. (|Bip shows that for a fixed C/b (or k), T is maxi- 



mal when E\\/Uh = 1 (see Fig. [TOl) . where E'j] can be 
tuned by changing the band edge V. 

In an FIF MTJ, for electrons in F with the same total 
energy E = e'I^^^^^ + ^T/i + ^'^'^ same transverse 
energy E±, the relative barrier height [Uh in Eq. (IBip ) 
is the same for spin up and spin down electrons. How- 
ever, the band edges (V|/|) are spin-dependent as in- 
dicated by the solid (spin up) and dashed (spin down) 
lines in Fig. [S) The longitudinal kinetic energies (i?|| in 
Eq. (|Bip ) is larger for spin-up than spin-down electrons. 
In our Stoner model we typically find /C/b ^ 1 whereas 

^ij/C/b > 2. According to Eq. (|BT]) and Fig. [TOl we find 
for the parameters in our Stoner model the surprising re- 
sult that T| < . In general, electrons close to the Fermi 
energy with small transverse wave vectors q show this in- 
verted polarization sign change behaviour. For large 
spin-up electrons tend to have higher transmission again, 
and the polarization becomes positive. 
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